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We report experimental results for convection near onset in a thin layer of a homeotropically 
aligned nematic liquid crystal heated from below as a function of the temperature difference AT 
and the applied vertical magnetic field H. When possible, these results are compared with theoretical 
calculations. The experiments were done with three cylindrical cells of aspect ratios (radius/ height) 
T = 10.6, 6.2, and 5.0 over the field range 8 < h = H/H F < 80 (H F = 20.9, 12.6, and 9.3 Gauss 
are the Freedericksz fields for the three cells). We used the Nusselt number AT (effective thermal 
conductivity) to determine the critical Rayleigh number R c and the nature of the transition. We 
analysed digital images of the flow patterns to study the dynamics and to determine the mean 
wavenumbers of the convecting states. For h less than a codimension-two field h ct ~ 46 the bifur- 
cation is subcritical and oscillatory, with travelling- and standing-wave transients. Beyond h c t the 
bifurcation is stationary and subcritical until a tricritical field h t = 57.2 is reached, beyond which it 
is supercritical. We analyzed the patterns to obtain the critical wavenumber k c and, for h < h c t, the 
Hopf frequency w c .In the subcritical range we used the early transients towards the finite-amplitude 
state for this purpose. The bifurcation sequence as a function of h found in the experiment confirms 
the qualitative aspects of the theoretical predictions. Even quantitatively the measurements of R c , 
k c and f2 c are reproduced surprisingly well considering the complexity of the system. However, the 
value of h c t is about 10% higher than the predicted value and the results for k c are systematically 
below the theory by about 2% at small h and by as much as 7% near h c t. At hct, k c is continuous 
within the experimental resolution whereas the theory indicates a 7% discontinuity. The theoretical 
tricritical field h\ h = 51 is somewhat below the experimental one. The fully developed flow above R c 
for h < hct has a very slow chaotic time dependence which is unrelated to the Hopf frequency. For 
hct < h < ht the subcritical stationary bifurcation also leads to a chaotic state. The chaotic states 
persist upon reducing the Rayleigh number below R c , i.e. the bifurcation is hysteretic. Above the 
tricritical field ht , we find a bifurcation to a time independent pattern which within our resolution is 
non- hysteretic. However, in this field range, there is a secondary hysteretic bifurcation which again 
leads to a chaotic state observable even slightly below R c . We discuss the behavior of the system in 
the high-field limit, and show that at the largest experimental field values R c and k c are within 6% 
and 1%, respectively, of their values for an infinite field. 



PACS numbers: 47.54.+r,47.20Bp,61.30-v 



I. INTRODUCTION 



Convection in a thin horizontal layer of an isotropic 
fluid heated from below by a heat current Qis well known 
as Rayleigh-Benard convection (RBC) When the 

fluid is a nematic liquid crystal (NLC), this phenomenon 
is altered in interesting ways. || NLC molecules are long, 
rod-like objects which are orientationally ordered rela- 
tive to their neighbors, but whose centers of mass have 
no positional order j||J. The axis parallel to the av- 
erage orientation is called the director h. By confining 
the NLC between two properly treated parallel plates 
H , one can obtain a sample with uniform planar (paral- 
lel to the surfaces, i.e. in the x-y plane) or homeotropic 
(perpendicular to the surfaces, or parallel to the z-axis) 
alignment of n. The alignment can be re-inforced by the 
application of a magnetic field H parallel to the intended 



direction of n. This is so because the diamagnetic sus- 
ceptibility is anisotropic, usually being larger in the di- 
rection parallel to the long axis of the molecules. The 
phenomena which occur near the onset of convection de- 
pend on the orientation of h and H . Q In this paper we 
are concerned with a horizontal layer of a homeotropi- 
cally aligned NLC in a vertical magnetic field (H = He z ) 
and heated from below. In that case, Q = Qe z is parallel 
to n when the system is in the conduction state. At a 
critical temperature difference AT C (H) the fluid will un- 
dergo a transition from conduction to convection. The 
precise value of AT C (H), the nature of the bifurcation at 
AT C (H), and the pattern- formation phenomena beyond 
AT C (H) are expected to depend in interesting ways upon 

A feature common to the homeotropic NLC and to an 
isotropic fluid heated from below is that the system is 
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isotropic in the horizontal plane. Thus the convection 
pattern may form with no preference being given to a 
particular horizontal axis unless the experimental appa- 
ratus introduces an asymmetry. In both cases, the con- 
vection is driven by the buoyancy force. However, the 
mechanism in the NLC case is more involved. The 
usual destabilization due to the thermally-induced den- 
sity gradient is opposed by the stiffness of the director 
field which is coupled to and distorted by any flow. Since 
relaxation times of the director field are much longer than 
thermal relaxation times, it is possible for director fluc- 
tuations and temperature/velocity fluctuations to be out 
of phase as they grow in amplitude. The existence of 
two very different time scales and this phase shift typ- 
ically lead to an oscillatory instability (also known as 
overstability) , i.e. the bifurcation at which these time- 
periodic perturbations acquire a positive growth rate is 
a Hopf bifurcation. |7^J[lJ This case is closely analogous 
to convection in binary-fluid mixtures with a negative 
separation ratio. Jl2|,[l3| In that case, concentration gra- 
dients oppose convection, and concentration diffusion has 
the slow and heat diffusion the fast time scale. As in the 
binary mixtures, the Hopf bifurcation in the NLC case 
is subcritical P,|l0||. For the NLC the fully developed 
nonlinear state no longer is time periodic. Instead, the 
statistically stationary state above the bifurcation is one 
of spatio-temporal chaos with a typical time scale which 
is about two orders of magnitude slower than the theo- 
retical inverse Hopf frequency. Jl4j However, it is possible 
to actually measure the Hopf frequency by looking at the 
growth or decay of small perturbations which are either 
deliberately introduced || or which occur spontaneously 
when the system is close to the conduction state and near 
the bifurcation point. 

A linear stability analysis of this system was carried 
out by several investigators . A very de- 

tailed analysis was presented by Feng, Decker, Pesch, 
and Kramer (FDPK) pcj . These authors also provided a 
weakly-nonlinear analysis, which allowed the distinction 
between sub- and supercritical bifurcations. Quantitative 
bifurcation diagrams were predicted for the nematic liq- 
uid crystal MBBA. In the present work we repeated and 
slightly extended the calculations for the material 5CB 
(see below) used in our experiments. Since the mate- 
rial parameters of MBBA and 5CB are similar, we found 
qualitatively the same bifurcation sequences as a func- 
tion of the field. Here we outline briefly the theoretical 
results and their relationship to our experimental results. 

As the magnetic field is increased, a subcritical Hopf 
bifurcation-line is expected to terminate at H ct in a 
codimcnsion-two point (CTP) beyond which the pertur- 
bations which first acquire a positive growth rate are 
at zero frequency. Close to but beyond the CTP this 
stationary bifurcation is predicted to also be subcritical. 
At an even higher field H t a tricritical point (TCP) is 
predicted to exist beyond which the primary bifurcation 
is expected to become supercritical. To our knowledge 
there are no predictions about the patterns which should 



form beyond either the subcritical Hopf bifurcation below 
the CTP or the subcritical stationary bifurcation between 
the CTP and the TCP. Although there are no explicit 
predictions of the patterns for H > H tl in analogy to 
isotropic fluids one might expect convection rolls above 
the supercritical bifurcation, unless non-Boussinesq ef- 
fects 1 18 19] yield a transcritical bifurcation to hexagons. 



The phenomena described above were previously ex- 
plored only partially by experiment. Except for recent 
measurements at relatively small fields [Q, the experi- 
ments have been qualitative or semi-quantitative |9|,|20|. 
In the present paper we report the results of an exten- 
sive sytematic experimental investigation of this system 
which covered a wide range of mag netic fields H . In 
agreement with previous work [p|,pO[, we find a subcrit- 
ical Hopf bifurcation at relatively small H which termi- 
nates in a codimension-two point (CTP). The CTP is 
located at a slightly higher field than the theoretical pre- 
diction. We measured the Hopf frequency lu c (H) from 
visualizations of the spontaneous small-amplitude early 
transients just above AT C . Except for the influence of the 
small shift of the CTP, we found ui c (H) to be in quanti- 
tative agreement with the theory. From these transients, 
we also determined the critical wavevector k c (H), and 
found it to be typically a few percent smaller than the 
theoretical value. The reason for this small difference 
between theory and experiment is not known. As re- 
ported previously [fl4|| , we found the convecting nonlin- 
ear state for AT above AT C to be one of spatio-temporal 
chaos (STC). Except at very small H, its characteristic 
wavenumber was smaller than k c and insensitive to H 
and AT. A long-time average of the structure factor of 
this state was consistent with the expected rotational in- 
variance of the system. Depending on H, the lower limit 
AT S at which this chaotic state made its hysteretic tran- 
sition back to the conduction state was found to be 10 
to 25% below R C (H). The convective heat transport was 
consistent with that of an isotropic fluid with an average 
conductivity given by X avg = (2Aj_ + A||)/3 where A_l and 
A || are the conductivities perpendicular and parallel to n 
respectively. This suggests a thorough randomization of 
the director orientations by the flow. 

Beyond the CTP we found a subcritical stationary 
bifurcation, as had been predicted |px| . The finite- 
amplitude state which evolved was also a state of STC, 
but beyond a certain field value greater than H ct it had 
distinctly different properties from the chaotic state at 
lower H. This difference was clearly evident from a 
discontinuity (as a function of H) of the characteristic 
wavenumber of the nonlinear state, which was larger at 
the larger fields. There is no theoretical guidance for the 
interpretation of these experimentally observed phenom- 
ena. 

The range H > H t was investigated for two cells with 
aspect ratios T — 6.15 and 5.01. We refer to them as 



cells 5 and 6 respectively (for details see Sect. Ill below). 
We found a primary bifurcation to a state with a hexag- 
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onal flow pattern. Within our resolution this bifurca- 
tion was non-hysteretic and the Nusselt number Af grew 
continuously from zero. The appearance of hexagons 
rather than rolls is attributable to non-Boussinesq ef- 
fects []l8|,[l9| which occur when the up-down symmetry 
is lifted by variations of the fluid properties over the cell 
height. Theoretically the bifurcation to hexagons is tran- 
scritical, and there should also be hysteresis associated 
with it. However, the hysteresis is often so small that 
it is unobservable even though hexagons are found over 
a substantial range. |2l| At sufficiently large AT/AT C 
the hexagons become unstable with respect to rolls. |19]] 
Since AT C decreases with the cell thickness d (~ d^*), 
the range of stability of hexagons should depend on the 
thickness of the fluid sample. However, for the thinner 
fluid layer of cell 5, the existence range of hexagons was 
limited by a different secondary instability and grew from 
zero very near the TCP to e = AT / AT C - 1 ~ 0.1 at the 
highest fields available to us. At this stability limit a 
hysteretic transition yielded the chaotic state, and sta- 
tionary rolls were never found. With decreasing AT, the 
chaotic state persisted down to AT S somewhat smaller 
than AT C . For the thicker fluid layer of cell 6, typically 
AT C was about 2 °C, and hexagons were found only up 
to e ~ 0.015 even at high fields. For e > 0.015, a pattern 
of rolls not exhibiting STC was observed, as expected for 
a weakly non-Boussinesq system. At even higher e and 
consistent with the measurements in cell 5, a hysteretic 
secondary bifurcation again yielded the chaotic state. 



II. PARAMETER DEFINITIONS AND VALUES 

The quantitative aspects of the instabilities are de- 
termined by four dimensionless parameters which are 
formed from combinations of the fluid properties [p2[ . 
They are ]lC| ] the Prandtl number 
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the ratio between the director-relaxation time and the 
heat-diffusion time 
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the Rayleigh number 
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and the dimensionless magnetic field 

h = H/ Hp 
with the Freedericksz field 
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In these equations is one of the viscosity coefficients, 
ac II is the thermal diffusivity parallel to n, k 33 is one of the 
elastic constants of the director field, \a is the anisotropy 
of the diamagnetic susceptibility, a is the isobaric ther- 
mal expension coefficient, and g is the gravitational ac- 
celeration. The time scale of transients and pattern dy- 
namics is measured in terms of the thermal diffusion time 



t v = d 2 /n\\ 



(6) 



Both h and R are easily varied in an experiment, and 
may be regarded as two independent control parameters. 
The availability of h in addition to R makes it possible to 
explore an entire line of instabilities. The parameters F, 
a, and t v are essentially fixed once a particular NLC and 
temperature range have been chosen, and even between 
different NLCs there is not a great range at our disposal. 
For 5CB at 25.6° (the material and mean temperature 
used in this work), we have a — 263 and F — 461. The 
value of t v is typically several minutes, but depends on 
the thickness of the fluid layer. It is given in the next 
Section for each of our cells. The critical value R c (h) of 
R and the fluid parameters determine the critical tem- 
perature difference AT C for a sample of a given thickness 
d. The realistic experimental requirement that AT C ~ 
a few °C dictates that the sample thickness should be a 
few mm. Typical values of Hp are 10 to 20 Gauss. Thus 
modest fields of a kGauss or so are adequate to explore 
the entire range of interest. 

In order to evaluate R c from AT C , h — H/Hp, and the 
theoretical values for R c (h),k c (h), and uj c (h), we used 
the material properties given in Ref. [p2[ . We followed 
closely the calculational methods of FDPK. In order to 
insure a sufficient resolution of any boundary layers we 
used Chebycheff modes in the Galerkin method (no more 
than 20 were required). 



III. EXPERIMENTAL APPARATUS AND 
SAMPLE PREPARATION 



The apparatus used in this work was described previ- 
ously |23|JT^ ]. We made measurements using three circu- 
lar cells of different thicknesses, identified as cells 4, 5, 
and 6. |24| The thickness and radius were d = 3.94 mm, 
r = 41.9 mm for cell 4, d — 6.60 mm, r = 40.6 mm for 
cell 5, and d — 8.88 mm, r = 44.5 mm for cell 6. The cor- 
responding radial aspect ratios T = r/d were 10.6, 6.15, 
and 5.01. The fluid was 4-n-pentyl-4'-cyanobiphenyl 
(5CB). All experiments were performed at a mean tem- 
perature of 25.6°C. The vertical thermal diffusion time 
was t v = 139, 383, and 694 s, and Hp was 20.1, 12.6, 
and 9.34 Gauss for cells 4, 5, and 6 respectively. De- 
spite the longer time scales involved for experiments in 
the thicker cells, cells 5 and 6 had an advantage over 
cell 4 due to the smaller field strengths and temperature 
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differences required to perform the measurements. To in- 
sure homeotropic alignment near the surfaces of the top 
and bottom plates of both cells, a surface treatment with 
lecithin was applied. 

Defect-free homeotropic samples were prepared by ap- 
plying a magnetic field while cooling the bath, and thus 
the sapphire top plate of the sample, from above the 
isotropic-nematic transition temperature Tjv/ to T < 
Tni- During this process, the bottom plate naturally 
lagged behind, and thus an adverse density gradient ex- 
isted. In the ncmatic/isotropic two-phase region even the 
relatively small thermal gradients associated with small 
cooling rates induced convection ]25|] . When the cool- 
ing was too rapid and the field too small this led to 
a nematic sample with defects which remained frozen. 
By using cooling rates of l°C/hour in the presence of a 
field of h > 17 over the temperature interval 36 to 34° C 
(Tjv/ = 35.1°C) and annealing at 34°C for an hour or two 
the defects healed and a defect-free homeotropic sample 
could be prepared. Further cooling could then be at least 
ten times as rapid without introducing new defects be- 
cause the threshold for convection in the nematic phase is 
large. Before each experimental run, the procedure was 
repeated. 

The critical temperature differences for the onset of 
convection were determined from heat-transport mea- 
surements. These are usually expressed in terms of the 
Nusselt number 



N = A e// /A|| 



(7) 



where An is the conductivity of the homeotropically 
aligned sample E3], and 



A e// = -Qd/AT 



(8) 



is the effective conductivity and contains contributions 
from diffusive heat conduction and from hydrodynamic- 
flow advection. Measurements of AT were made by de- 
termining the heat current Q required to hold AT con- 
stant. At each AT, the heat current and temperature of 
the bath and bottom plate were measured at one-minute 
intervals for three to five hours, when typically all tran- 
sients had died out. 

In addition to heat-flow measurements, we also visual- 
ized the convective flow patterns. The homeotropic sam- 
ples were translucent even for d as large as several mm. 
It was just about possible to see features of the bottom 
plate in typical ambient lighting. Any director distortion 
by convection rolls or domain walls generated opaque re- 
gions with enhanced diffuse scattering which were easily 
visible. It should be kept in mind that the optical sig- 
nal in the images has a complicated relationship to the 
hydrodynamic flow fields, and that quantitative infor- 
mation about velocity- or temperature-field amplitudes 
could not be obtained. Such quantities as the wavevector 
of the patterns or frequencies of traveling convection-rolls 
could of course be determined quantitatively. 



The samples were illuminated from above by a circular 
fluorescent light. Digital images were taken from above 
by a video camera which was interfaced to a computer. 
Typically 50 to 200 images were averaged to improve the 
signal-to-noise ratio. Averaged images were divided by 
an appropriate reference image to reduce the influence 
of lateral variations in illumination and of other optical 
imperfections. Some images were processed further by 
filtering in Fourier space. 
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FIG. 1. Examples of the Nusselt number as a function of 
AT for cell 5. The upper figure is for h — 15 and the lower one 
is for h = 50. Open circles were taken with increasing and 
solid circles with decreasing AT. The transitions between 
conduction and convection are indicated by the arrows. 



IV. RESULTS 

A. Nusselt Numbers and Critical Rayleigh Numbers 

Figure [l] shows J\f for cell 5 as a function of AT for 
two field strengths h = 15 and h — 50. The open cir- 
cles were obtained with increasing, and the solid circles 
with decreasing AT. For the lower fields (h < 20) N 
decreased below one when convection started. This can 
be understood because the convecting sample has a dis- 
torted director with a component perpendicular to Q. 
The contribution from this component to the conductiv- 
ity corresponds to Aj_, which is less than the conductiv- 
ity An of the homeotropic case. |^3| It turns out that for 
small fields the direct hydrodynamic contribution to the 
heat flux is smaller than the decrease in the heat flux 
due to the director distortion by the flow. For the higher 
fields (h > 35), N remained above one in the convect- 
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ing state. Thus with the higher fields the hydrodynamic 
contribution to the heat flux is greater than the decrease 
in the heat flux due to any distortion of the director. 
Both examples in Fig. Ill demonstrate the predicted and 
previously observed |]jlJ,|2(J hysteretic nature of the bi- 
furcation, i.e. as AT was decreased, the conduction state 
was reached at a value of AT equal to AT S < AT C . 

From data like those in Fig. |l|, critical temperature dif- 
ferences AT C were determined with an uncertainty of less 
than one percent. The corresponding Rayleigh numbers 
are shown in Fig. | as a function of h 2 . The open cir- 
cles were obtained in cell 4, the filled ones in cell 5. The 
good agreement between the two data sets confirms the 
expected scaling of the field with Hp. It also shows that 
using the fluid properties at the mean temperature does 
not lead to systematic errors in R c even for cell 4 where 
AT C is over 10°C. One sees that R c is quadratic in h at 
small h, as is expected because the system should be in- 
variant under a change of the field direction. The solid 
line follows the theoretical prediction and agreement be- 
tween theory and experiment is excellent. 
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FIG. 2. Critical Rayleigh numbers for the onset of convec- 
tion as a function of h 2 . Open and filled symbols were ob- 
tained in cells 4 and 5, respectively. The line is the theoretical 
prediction. 

Results for R c over our full experimental field range are 
shown as a function of h in Fig. ^. Here we include data 
taken with cell 6 as open squares. The data for R c reveal 
a sharp maximum at h = 44.3. We interpret this field 
value as the codimension-two point h ct and indicate it in 
Fig. H by the dashed vertical line. The solid line in the fig- 
ure is the theoretical prediction for R c , evaluated for the 
properties of our sample. For the entire range h < h ct , 
the theoretical result is in excellent agreement with the 
data. However, the theory gives h c t = 41.8 which is 
slightly lower than the experimental value. Above h c t 
the measurements of R c are systematically larger than 
the calculation, although the largest discrepancy is only 
about 4%. 

The triangles in Fig. [|show the lower limit of existence 
(the "saddle- node" Rayleigh number R s ) of the finite- 
amplitude convecting state as determined from data like 
those in Fig. [|. They suggest that the tricritical bifur- 



cation is located near h = 59, which is larger than the 
theoretically calculated value h t ~ 51. However, we will 
return later to the best estimate of ht- 

Measurements similar to those shown in Fig. |§| were 
made by Salan and Fernandez- Vela (SF), using the 
nematic liquid crystal MBBA. Their results are shown in 
Fig. |], together with the theoretical curve for that case 
pq | . The data and the curve illustrate that there are sig- 
nificant quantitative differences between the bifurcation 
lines of different nematics. In Fig. || the experimental 
points lie on average about 25% above the theoretical 
curve, and the lower hysteresis limit is further below the 
bifurcation line than we found for 5CB. 
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FIG. 3. Critical Rayleigh numbers R c and saddle-node 
Rayleigh numbers R s over the experimentally accessible field 
range. The open circles, filled circles, and open squares are 
R c for cells 4, 5, and 6 respectively. The open and filled tri- 
angles are R s for cells 4 and 5 respectively. The dashed line 
indicates the location of the codimension-two point as found 
experimentally. The solid line is the theoretical prediction for 
R c - The plusses at h — 35 were obtained with short equili- 
bration times and cell 5 (see text). 
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FIG. 4. Results for R c and R a from Ref. Q obtained with 
MBBA. The theoretical curve for R c was computed from typ- 
ical fluid properties of MBBA @ . 

The equilibration times after each temperature step used 
by SF were 30 minutes, which is a factor of six to ten 
shorter than those of our experiments. In addition the 
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temperature steps of SF were a factor of two larger than 
ours, yielding a difference in the average rate of change 
of the temperature of a factor of 12 or more. Looking for 
an explanation of the difference between the experimen- 
tal and theoretical R c revealed in Fig. [|, we conducted 
one run with equilibration times similar to those of SF, 
but using our 5CB sample. It gave the plusses in Fig. [|. 
As can be seen, these results do not differ significantly 
from the data taken with our usual longer equilibration 
times. Thus we have no explanation for the difference be- 
tween the SF data for R c and the theoretical curve. How- 
ever, the agreement between our runs with the different 
equilibration times implies that our usual experimental 
procedure yielded quasi-static results. 




FIG. 5. A sequence of images of the travelling/standing 
wave transients for h — 32 and e = 0.015 in cell 5. The num- 
ber in each image corresponds to the time, in units of t v , that 
elapsed since flow first became visible. A time series of the 
pixel intensity was taken at the location marked in the top 
left image. 

B. Hopf Frequency and Critical Wavevector 

Once the critical Rayleigh numbers were measured, a 
detailed analysis of the patterns could be undertaken. 
At first we will characterize the Hopf bifurcation for h 
below h c t- Since in that field range the bifurcation is 
subcritical, we had to use the small-amplitude transients 
to determine to c and k c . Figure ^ shows images from 
cell 5 which are characteristic of these patterns. They 
were taken at the times (in units of t v = 383 s) indi- 
cated in each figure after the pattern initially became 
visible. This typically occurred around one hour after e 
was raised from below zero to around 0.015. Inspection 
of successive images revealed that the transients could be 
either traveling or standing waves, sometimes with both 



occurring at different locations in the same cell. In the 
top left image of Fig. ||, a location is indicated at which 
a time series of the pixel intensity was acquired. This 
time series is shown in Fig. |^ along with its correspond- 
ing power spectrum. The length of the time series was 
limited by the rapid growth of the pattern to its finite- 
amplitude steady state. Because of this, only a small 
number of periods could be obtained before the finite- 
amplitude state was reached. Thus to avoid errors asso- 
ciated with incommensurate sampling, the data was win- 
dowed before its Fourier transform was evaluated. This 
process was repeated at several pixel locations in the cell. 
The signal from the second harmonic was often found to 
be stronger than that from the fundamental. Thus it was 
used to calculate the frequency. The frequencies at dif- 
ferent locations generally were within a percent of each 
other, and were averaged to determine the critical Hopf 
frequency lu c . 
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FIG. 6. The upper figure is the time series of the pixel in- 
tensity at the location shown in the top left image of Fig.^. 
The lower figure is the power spectrum of that time series. 
The mean frequency calculated from the solid circles was used 
to determine the Hopf frequency. 



The dependence upon h of the measured w c is com- 
pared with theory in Fig. [?]. The arrow indicates the lo- 
cation of the theoretical codimension-two point while the 
dashed line represents the experimental determination of 
h ct . As can be seen, away from the codimension-two 
point the agreement with the measurements is excellent. 
In accordance with theory, the experimental uj c changes 
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discontinuously to zero at h ct , above which the bifurca- 
tion is stationary. 

By evaluating the Fourier transforms of images such as 
those in Fig. |^, the critical wavenumber k c of the patterns 
could be measured. The transforms were based on the 
central parts of the images by using the filter function 
W(r) = {1 + cos[(n)(r/r )]}/2 for r < r and W(r) = 
for r > r Q . Here r a was set equal to 85% of the sample 
radius. Time averaging the square of the modulus of the 
transforms over the length of the time series yielded the 
structure factor < S(k) >. Figure || shows the azimuthal 
average < S(k) > of < S(k) > for the run at h = 32. We 
used a weighted average of the three points nearest the 
peak of the second harmonic of < S(k) > to calculate k c . 
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FIG. 7. The Hopf frequency ui c as a function of h. Open 
and filled circles are for cells 4 and 5, respectively. Open 
squares are for cell 6. The solid line is the theoretical predic- 
tion for LOc- The dashed line (arrow) indicates the location of 
the codimension-two point as found experimentally (theoret- 
ically) . 




FIG. 8. The azimuthal average of the time-averaged struc- 
ture factor < S(k) > of the travelling/standing wave tran- 
sients at h = 32 for cell 5. The mean wavenumber calculated 
from the solid circles was used as the critical wavenumber. 

Figure || displays the results for k c for all h together 
with the theoretical analysis. For h < h ct the measured 
critical wavenumber of the transients is systematically 
smaller than the theoretical one. When the codimension- 



two point is approached, the experimental wavenumbers 
make a smooth rather than discontinuous transition to 
those associated with the stationary bifurcation, whereas 
the theory predicts a 7% discontinuity of k c at h ct . The 
reason for these discrepancies is as yet unknown. Above 
h ct , the agreement between the experimental and theo- 
retical wavenumbers is excellent. 

As shown explicitly for R c in Fig. R c ,uj c , and k c are 
proportional to h 2 for small h. 
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FIG. 9. The characteristic wavenumbers of the observed 
patterns as a function of h. Circles: The wavenumber of the 
small-amplitude transients. Triangles: The wavenumber k n 
of the fully developed spatially and temporally chaotic flow, 
as measured close to the onset of such flows. Open and filled 
symbols were obtained in cells 4 and 5, respectively. 

C. Nonlinear State below the tricritical field h ct 

Because of the subcritical nature of the bifurcation for 
h < h ct a finite-amplitude state develops directly at on- 
set. The time dependence and spatial structure of this 
state are very different from that of the small-amplitude 
transient state. The first two rows of Fig. show typical 
images of the patterns from cell 5 which are characteristic 
of the fully developed flow. They are from a single ex- 
perimental run with constant external conditions. They 
were taken at the times indicated in each image, in units 
of t v — 383 s, which had elapsed since e had been raised 
from below zero to 0.014. The convection rolls of the 
fully developed flow show an irregular time dependence, 
with typical time scales around a hundred times longer 
than the inverse Hopf frequencies of the transients. One 
can see that the "chaotic" behavior is associated with the 
formation of defects and the continuous reorientation of 
the convection rolls. This continuous reorientation of the 
rolls is evident in the rightmost image in the bottom row 
of Fig. |l^ (labeled "Avg"). It shows the time average 
of the structure factor. The average involved 75 images 
taken over a total time period of 724t v (over three days) . 
It is seen to contain contributions at all angles, consis- 
tent with the idea of a statistically stationary process of 
non-periodic pattern evolution and with the expected ro- 
tational symmetry of the system. Similar results for cell 
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4 have been shown previously, uA 




FIG. 10. Top two rows: a sequence of images taken with 
constant external conditions (h — 50, e = 0.014) for cell 5. 
The time elapsed since e was raised from below zero (in units 
of t v = 383s) is given in the top left corner of each image. 
Bottom row: The structure factor of two of the images shown 
above, and the average of the structure factor of 75 images 
spanning a time interval of 724t„. The structure factor was 
obtained using a Hanning window, and thus is dominated by 
the patterns near the cell center. 



FIG. 11. A temporal succession of images during the tran- 
sient leading from conduction to convection when AT was 
raised slightly above AT C for cell 5. The field was h = 50. 
The numbers are the elapsed time, in units of t v , since the 
threshold was exceeded. 



When h was increased above h c t, the nature of the pat- 
tern at hrst did not change noticeably. For instance, as 
evident in Fig. ^, the characteristic wavenumber of the 
pattern (as denoted by the triangles) remained close to 
3.4 for h < 55. Over this field range the patterns of the 



fully developed flow look similar to those illustrated in 
Fig. HQ, i.e. they exhibit spatio-temporal chaos. 



CD 

E 



1.4 
1.2 
1 

1.4 
1.2 
1 

1.4 
1.2 
1 

1.4 
1.2 
1 

1.4 
1.2 
1 

1.4 
1.2 
1 





h = 55 



-o- oq - 



-oeGr - 



35 



I.. 



■ eeeor - 



25 



3-O-O-O-Ot 



15 



- e-o-e -°j<5^b»c*e ,0 * : 



-0.3 



-0.2 



-0.1 

AT / AT C - 1 



0.1 



FIG. 12. Nusselt-number measurements for cell 5 illustrat- 
ing the variation of the size of the hysteresis loop between 
conduction and convection with h. The number in the upper 
left corner of each plot is the held h. Open circles were taken 
with increasing and solid circles with decreasing AT. The 
arrows show the values of e a = R 3 /R c — 1. 



It is instructive to examine the transients which 
lead from the small-amplitude to the finite-amplitude 
statistically-stationary state. This is done in Fig. O. 
Here the number in each image gives the time, in units 
of t v , which had elapsed since AT was raised slightly 
(1.4%) above AT C . At t = 4.7 small-amplitude tran- 
sients like those in Fig. || are evident in part of the 
cell. By t = 9.4 these had filled the cell and grown 
to a saturated amplitude. At this stage they formed 
nearly-straight parallel rolls with a wavenumber which 
was smaller than k c . However, these straight rolls turned 
out to be unstable to a zig-zag instability. In the end, 
this instability led to the spatially and temporally disor- 
dered pattern as shown in Fig. [10| Thus, we see that a 
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secondary instability led to a chaotic state rather than to 
a new time-independent pattern. This phenomenon most 
likely is similar to the one encountered in very early ex- 
periments on spatio-temporal chaos using liquid helium 
p7|p8[ , where ordinary RB convection became chaoti- 
cally time dependent, most likely because the secondary 
skewed- varicose instability p9| was crossed. 
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FIG. 13. Nusselt-number measurements for cell 5, normal- 
ized by the averaged conductivity \ avg = (2Ax + Am)/3, 
as a function of the Rayleigh number Ravg computed with 
Kavg = ^avg/pC'p. The data are for h = 15 (open circles), 
25 (solid circles), 35 (open squares), 45 (solid squares), 50 
(open triangles), and 55 (solid triangles). The diamonds are 
for h = 64, which is in the supercritical region. Here the solid 
diamonds are for the hexagons/rolls which form supercriti- 
cally, and the open ones are for the chaotic finite-amplitude 
state which forms via a secondary bifurcation (see Fig. H be- 
low). The cross corresponds to the critical Rayleigh number 
R c = 1708 of an isotropic fluid. 



Heat-transport measurements of the fully developed 
flow are shown in Fig. |l2| for several field values as a 
function of e. They illustrate the evolution with h of the 
hysteretic nature of the bifurcation. As can be seen also 
in Fig. ||, the hysteresis |e a | increased with h for h < h ct 
from about 10% at the low fields to nearly 25% close to 
the codimension-two point. Above this point the hystere- 
sis decreased and suggested the existence of a tricritical 
point near h t ~ 59 (see below for more detail about the 
tricritical region). 

When the Nusselt numbers in Fig. [l^ are plotted 
against the Rayleigh number R, they fall on or approach 
a single curve independent of h. This suggests that the 
convection in the chaotic state is sufficiently vigorous to 
achieve nearly complete randomization of the director 
orientations, regardless of h. In that case one would 
expect that the system should behave approximately 
like an isotropic fluid, with an averaged conductivity 
Aat) 9 = (2Aj_ + A||)/3. Thus we plot in Fig. [ll| a modi- 



fied Nusselt number Af given by the ratio of the effective 
conductivity of the convecting state to X avg as a function 
of Ravg, where R avg is computed using K avg = Kvg/pCp 
in Eq. || rather than km. At all but the highest fields 
(where the primary bifurcation is supercritical) the data 
reach the common curve. At small R a vgi this curve ex- 
trapolates to J\f = 1 near R avg = 1708 (the cross in 
the figure), which is the critical Rayleigh number of an 
isotropic fluid. An analogous behavior has been observed 
in binary-mixture convection with negative separation 
ratios 'F, 1 30 1 where the bifurcation is also subcritical. In 
that case the convective flow achieves thorough mixing of 
the concentration field and Af approaches a curve which 
is independent of In both cases the mixing achieved 
by the flow can persists because of the existence of a 
slow time scale, namely that of director or concentration 
relaxation. 

Further support for the idea that the chaotic flow in 
some respects can be approximated by isotropic-fluid 
convection is found in Fig. ^, where for h ^$ 55 the 
wavevectors (triangles) are independent of h and much 
closer to the critical value k l c so = 3.117 than to the critical 
values k c (h) of the anisotropic system (circles in Fig. ||). 
Exact agreement with k l c so would of course not be ex- 
pected even for a genuine isotropic fluid because of the 
finite flow amplitude and various wavenumber-selection 
processes. 

Lastly we note that an extrapolation of the data in 
Fig. [l3| to Af = 1 and R avg = 1708 yields an initial slope 
Si of Af - 1 = Si(R avg /1708 - 1) of about 0.6. For a 
laterally infinite system of straight rolls in an isotropic 
fluid with a large Prandtl number one expects S\ ~ 1.43. 
ifjlf However, experiments in finite systems with mod- 
est aspect ratios |}2| have always yielded smaller values, 
usually in the range of 0.6 to 1. Particularly when many 
defects are present, as in our case, one would expect the 
heat transport to be suppressed relative to that of a per- 
fect straight-roll structure. 



D. Tricritical Region and beyond 

This section is devoted to the phenomena which oc- 
cur near the tricritical field h — h t . At first the Nusselt 
numbers and the patterns are described and the corre- 
sponding bifurcation diagram is given. Further subsec- 
tions deal with the precise determination of ht and with 
hexagons observed near threshold for h > h t . 



1. Nusselt Numbers and Patterns 

From the measurements of the Nusselt number (see 
Fig. [12]) there is clear evidence of a tricritical field ht, 
above which the primary bifurcation is supercritical. For 
instance for h ~ 60, measurements of Af revealed no hys- 
teresis at the primary bifurcation and within our resolu- 







tion Af grew continuously from one beyond AT C . This is 
exemplified for cell 5 and h — 63 in Fig. The open 
(solid) circles correspond to the stable states reached 
by increasing (decreasing) AT. |}3| This behavior of TV 
stands in sharp contrast to that shown in Fig. ^| for lower 
fields. 




FIG. 14. Nusselt-number measurements for h = 63 and 
cell 5 illustrating the supercritical nature of the bifurcation 
characteristic of the high fields. Open circles were taken with 
increasing and solid circles with decreasing AT. 
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FIG. 15. Bifurcation diagram in the region of the R — h 
plane close to the tricritical point. Solid circles: primary 
bifurcation. Solid triangles: R s . Solid squares: hysteretic 
secondary bifurcation to chaotic convection. Shaded areas la- 
beled k n — 3.4, 3.6, and 3.9 correspond to chaotic regimes 
with different mean wavenumbers. The wedge-shaped area 
labeled Rolls/Hex shows the parameter range over which 
time-independent convection is stable. For cell 5, the pat- 
tern is hexagonal in this entire region. For cell 6, the pattern 
is hexagonal in this region for e < 0.015. For larger e but still 
in this region it consists of time-independent rolls. 



Besides the Nusselt number the analysis of the patterns 
gives important additional insight in particular with re- 
spect to secondary transitions. In Fig. [15] all the available 
information has been condensed in a bifurcation diagram 
for the vicinity of the tricritical point. 




FIG. 16. Image of the hexagonal flow in cell 5 for h = 65 
and e = 0.01. The pattern was essentially the same over the 
entire existence range of hexagons. 



At first we will focus on the wedge labeled Rolls/Hex 
where time-independent convection is stable. For both 
cell 5 and 6, a seemingly supercritical primary bifurca- 
tion led to a hexagonal pattern. For cell 5 this pattern is 
shown in Fig. [ll| The range of e over which the hexagons 
were stable differed in the two cells. In cell 5, he xag ons 
remained stable up to R n (h) (solid squares in Fig. [15]) for 
the entire range of h. At R n (h) a transition to a spatially 
and temporally chaotic roll pattern with a lower charac- 
teristic wavenumber occurred. For h < 75 this transition 
was distinguished by a jump in Af as well as the onset of 
time dependence of A/", as illustrated in Fig. [lj]. The up- 
per figure gives the steady-state Af and shows the jump 
at e„ = R n /R c — 1. The lower figure is the time series of 
Af obtained in the same run. Here e was held constant 
for a five-hour period at each of the eleven successively 
increasing values. The data show that Af is steady below 
and time dependent above e n (h). For ft, > 75 the discon- 
tinuity in Af was no longer pronounced, but a transition 
to time dependence still occurred at e n (h). Thus, de- 
pending on the field, either of these two indicators was 
used to determine the location of e n (h). In the thicker 
cell 6, a transition from hexagons to rolls occurred near 
e ~ 0.015, independent of field strength. This transi- 
tion was not associated with a measurable change in the 
wavenumber, a jump in Af, or a time dependence of Af. 
A further increase of e again led to a transition at e n {h) 
from steady rolls to the chaotic state, consistent with 
the cell 6 experiments. The results for e n (h) obtained in 
cells 5 and 6 are shown in Fig. [l^ as circles and triangles 
respectively. 

On the basis of the usual Landau equation for a tricrit- 
ical bifurcation, one would expect the hysteresis to grow 
gradually as h is reduced below h t . However, within our 
resolution this was not the case and a hysteretic primary 
bifurcation to a chaotic state occurred immediately be- 
low ht- Indeed, the secondary bifurcation line e n (h) for 
h > ft t met the primary bifurcation line at h t within ex- 
perimental resolution, as can be seen already in Fig. H5] 
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FIG. 17. Nusselt-number measurements for h — 61 in cell 
5. The upper figure illustrates the dependence of the Nus- 
selt number on e. The dotted line indicates e„, where the 
transition from hexagons to rolls occurred when AT was in- 
creased. The lower figure illustrates the dependence of the 
Nusselt number on time. Time is measured in units of 5 
hours, i.e. the time between steps in e. 
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agree with the tricritical field obtained from the slope 
of the Nusselt number (see the next Section). The bifur- 
cation for h > ht is supercritical but the amplitude at 
constant e > diverges as h approaches h t from above. 
Since secondary bifurcations occur at finite values of the 
amplitude, we expect e n (h) to vanish at h t . Therefore 
the e n measurements provide a relatively precise lower 
limit ht = 57.6 for the tricritical field. 
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FIG. 19. Nusselt-number measurements for h — 64 in cell 
5. Open circles were taken with increasing and solid circles 
with decreasing AT. The dotted line indicates e n . 

The secondary bifurcation at e n is strongly hysteretic. 
As illustrated in figure |l9|, when e was decreased from 
e > £n, a transition back to hexagons did not occur. In- 
stead, the chaotic state persisted to values of e slightly be- 
low zero, at which point the conduction state was reached 
(see also the solid triangles in Fig. 15). In a separate sec- 
tion we will come back to the hexagons. 
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FIG. 18. Values of e n (h) where a transition to a spatially 
and temporally chaotic state occurred. Circles and triangles 
were obtained in cells 5 and 6, respectively. The lines repre- 
sent a fit of a quadratic polynomial in (h/h n — 1) to the data. 
The fits extrapolate to zero at h n = 58.3 ± 0.7 for cell 5 and 
h n = 58 ± 3 for cell 6. 

It is shown more explicitly in figure |l8|, where the range 
e n of time-independent patterns vanishes near h = 58. 
Fitting a quadratic polynomial in (h/h n — 1) to e n (h) 
yields h n = 58.3 ±0.7 (cell 5) and h n = 58 ±3 (cell 6) for 
the field where e n vanishes. Within error, these values 



2. Determination of the Tricritical Point 

As h approaches the tricritical point from above, the 
initial slope S\,r of J\f for rolls is expected to diverge as 
l/(h— h t ). For cell 6, we estimated Sx >r (h) from data for 
e > 0.015 where rolls were observed. At a given h, Si. r 
was determined by fitting the polynomial 

N= l + S 1 , r e + S 2 , r e 2 (9) 

with e = AT/AT C - 1 to the data. The parameters AT C , 
Si jr , and i!?2,r were adjusted in the fit. Figure shows 
the dependence of 1/Si, r upon h as solid circles. The 
fitting procedure did not yield highly accurate values be- 
cause the Nusselt data for e < 0.015 had to be excluded; 
thus the error bars for 5i jr are relatively large. The line 
is a fit of a quadratic polynomial in (h/h t — 1) to the re- 
sults for l/Si }r (h). This fit indicates the tricritical point 
to be at ht =57.2 ±2.6. 

In order to compare the theory with the measurements 
above ht, we calculated Si <r using the properties of 5CB. 
The results for 1/Si jr are shown as a dashed line in 
Fig. [2(]. There is quite reasonable agreement with the 
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experimental data for h > 60, particularly when it is 
considered that the initial slope of M in finite systems 
usually is smaller than the theoretical value for the infi- 
nite system. However, the theory yields a tricritical field 
h\ h — hi which differs significantly from the experimen- 
tal estimates. Wc note that this difference is in the same 
direction as and somewhat larger than the corresponding 
one for the codimension-two point. We have no explana- 
tion for this difference. 
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FIG. 20. The dependence on h of the reciprocal of the ini- 
tial slope Si of the Nusselt number as obtained from a fit of 
the data to Eq. ^| The open circles are 1/Si,h for cell 5, and 
the filled ones are 1/Si, r for cell 6. For cell 5, data over the 
range < e < e„ were used, and the pattern was hexagonal. 
For cell 6, data over the range 0.015 < e < e n were used, and 
the pattern was one of rolls. The solid line represents a fit of a 
quadratic polynomial in (h/ht — 1) to the cell 6 5i, r data. The 
field ht = 57.2 ± 2.6, where for cell 6 extrapolates to 

zero, is interpreted to be the tricritical point. The theoretical 
results for l/5i, r are given by the dashed line. They yield a 
tricritical point at h\ h = 51. The theoretical results for 1/Si,h 
are given by the dash-dotted line. The location of the codi- 
mension-two point is given by the solid (experimental) and 
dashed (theoretical) short vertical lines. 



3. Hexagons 

In this section we will discuss in more detail the hexag- 
onal patterns. Hexagonal patterns at onset may be at- 
tributable to departures of the physical system from the 
Oberbeck-Boussinesq (OB) approximation p8lff9|, i.e. to 
a variations of the fluid properties over the imposed tem- 
perature range. For isotropic fluids it has been shown 
that non-OB effects lead to hexagons at a transcritical 
(hysteretic) primary bifurcation fl9] , |2~ij| . Below onset, for 
e a < e < 0, both hexagons and the conduction state are 
stable. Above onset, hexagons are stable for < e < e r . 
For e r < e < €b hexagons and rolls are both stable, while 
for e > €\, only rolls are stable. When the thickness of the 
fluid layer is increased, AT C is reduced and thus depar- 
tures from the OB approximation become smaller. Thus 
the range of e over which hexagons are stable is reduced 
when the thickness of the fluid layer is increased, as seen 



in the experiment by comparing cells 5 and 6. 

A stability analysis of RBC with non-OB effects in a 
homeotropically aligned NLC has not yet been carried 
out and would be very tedious. Thus, in order to obtain 
at least a qualitative idea of the expected range of stable 
hexagons, we used the theoretical results for the isotropic 
fluid with the fluid properties of 5CB. The values of e a , 
etc. are determined by a parameter V which was de- 
fined by Busse []l9| and is given by V = £|_ 7,-'Pj , with 
7o = -Ap/p , 7i = Aa/2a , 72 = Avjv , 73 = AA/A , 
and 74 = ACp/Cp . Here p is the density, a the ther- 
mal expansion coefficient, v the kinematic viscosity, A the 
conductivity, and Cp the heat capacity. The quantities 
Ap, etc. are the differences in the values of the properties 
at the bottom (hot) and top (cold) end of the cell. For 
A we used Ay, and for v we used ct^jlp. The coefficients 
Vi in the equation for V are given by Busse. [jl9| How- 
ever, here we use the more recent results p4| for large 
Prandtl numbers P = 2.676, Vi = -6.631, V% = 2.765, 
V 3 = 9.540, and V 4 = -6.225 where V 3 differs signifi- 
cantly from the earlier calculation. 

At the fields where the hexagons were observed, the 
temperature difference across cell 5 (cell 6) was close to 
5.05°C (2.07°C). At these temperature differences, we 
obtained V = -1.5, e a = -1.7 x 10~ 4 , e r = 1.5 x 10~ 2 , 
and e;, = 5.3 x 10~ 2 for cell 5. For cell 6 the values 
are V = -0.6, e a = -2.8 x 10~ 5 , e r = 2.5 x 10~ 3 , 
and eb — 8.7 x 10~ 3 . From these estimates, it follows 
that the hysteresis of size e a is too small to be notice- 
able in either cell with our resolution. The largest value 
of e at which hexagons could exist in cell 6 would be 
e = 8.7 x 10~ 3 . However, we observe hexagons to ex- 
ist to nearly twice this value. If the same is true for 
the thinner cell, the hexagon-roll transition attributable 
to non-OB effects should happen at a value of e greater 
than e n (h) for the field range over which the experiments 
were performed. Thus, instead of leading to a time in- 
dependent state, as observed in cell 6, the hexagon-roll 
transition is preceded by a transition to a state exhibiting 
spatio-temporal chaos at e n . 

The hexagonal pattern may be regarded as a superpo- 
sition of three sets of rolls with amplitudes A t , i = 1,2,3, 
corresponding to the three basis vectors at angles of 120° 
to each other. According to the Landau model [|l3| the 
steady-state amplitudes Ai are determined by 

eAt + h(A 2 2 + A3) - gil A\ - g 12 A 1 Al - g^Aj = 

(10) 

and the corresponding cyclic permutation for i = 2,3. 
Since all amplitudes are expected to be equal in hexagons 
{Ax — A 2 = A 3 = A), one has 

eA + bA 2 - (g + 2g)A 3 = (11) 

where g is the self-coupling coefficient gu and g is the 
cross-coupling coefficient g\2 = 313. Because of the term 
bA 2 , the bifurcation is transcritical and thus hysteretic. 
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However, as we discussed above and as is shown for in- 
stance by the data in Fig. [l4[ this effect is not resolved 
in the experiment because the coefficient b (which is de- 
termined by V) is too small. Thus we neglect the term 
bA 2 , and have to a good approximation 



N - 1 = 3 A 2 



3c 



2ff 



(12) 



At the tricritical point g vanishes as g = ga{h — h t ). 
However, there is no reason why g should vanish also at 
h t . Thus, one would expect the slope Si t h — 3/(<? + 2g) 
of Af near e = to remain finite at h t and equal to 3/2g. 
To test this idea, we fitted Af for cell 5 over the e range 
where hexagons were observed (i.e. up to e n ) to an equa- 
tion like Eq. |[ Since data quite close to threshold could 
be used, the results for Si t h are much more precise than 
those for cell 6. They are given in Fig. ^0|as open circles. 
One can see that 1/ S\ t h is non-zero at h t . It extrapolates 
to zero near h = 53, which is well below ht = 57.2 ± 2.6. 
Unfortunately the relatively large uncertainty of ht and 
Si^ r prevents the accurate determination of g. At ht, we 
find g = 3/2Si t h — 0-3. With increasing h, g also in- 
creases. For instance, the data in Fig. suggest that 
g = 3/2S hh - l/25i, r ~ 0.8 for h = 66. 

The experimental results for 1/ S\j L cannot agree quan- 
titatively with the theory because we already know that 
hf 1 is lower than the experimental value. Nevertheless 
we calculated 1/Si t h, and give it in Fig. as the dash- 
dotted line. We see that the relationship between 1/Si.h 
and l/Si >r is quite similar in theory and experiment. 



E. The High-Field Limit of R c and k c 

It is highly probable that there exists a high-field 
regime where the convection phenomena become inde- 
pendent of the field since the director is then frozen in 
the homeotropic configuration. It is instructive to exam- 
ine whether the experimental data extend to sufficiently 
high fields to fully reveal this behavior. Building on the 
results of FDPK pToj ] , one can show that in a one-mode 
approximation the neutral curve in the limit h — > oo is 
given by 



R (k) 



-hk 2 



(A^/A||)fc 2 + ^ 2 

2k 2 1 2 
2r]i 2t]2 2ai 

CI4 Q!4 Qf4 



Qf4 



(13) 



Here rji = (a 4 + a 5 - a 2 )/2 and 772 = (013 + «4 + 
ao}/2 are Miesowicz viscosity coefficients, g and I\ = 
0.69738, 1 2 = -12.3026, and b = 1.5056tt. This leads 
to R c ,oo(k c ) = 3090.6 and fc Cj00 = 4.294. Using many 
modes, one obtains numerically R c ,ao{k c ) = 3056.6 and 
k c , 00 = 4.328, close to the one-mode result. One can see 
that only the viscosities enter into i? Cj00 , and not the elas- 
tic constants. This is so because the director is held rigid 



by the field. R c .oa is larger than the isotropic- fluid value 
because of the additional viscous interaction between the 
flow and the rigid director field. 

In the high- field limit, we obtain 



R c — R c 



Rx/h 2 + 0{h- 



(14) 



The coefficient Ri has not been calculated in detail, but 
is proportional to fejjjf . The fact that at order 1/h 2 elastic 
constants enter suggests the beginning of some director 
distortion by the flow. 
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FIG. 21. The critical Rayleigh number as a function of 
1/h 2 . The solid line was calculated usine the parameters of 
Ref. The open triangle is R c ,oc from Eq. |l3[ The solid 

triangle is the many-mode numerical result for R c ,oc- The 
solid circles are from cell 5, and the open squares are from 
cell 6. The vertical bars indicate the location of the codi- 
mension-two point (right bar) and the tricritical point (left 
bar). The plusses, crosses, and dashed line show what hap- 
pens to the data and the theoretical curve if the viscosity «4 
is increased by 7.5%. 

In Figs. |2l] and [2^ we show the experimental data and 
theoretical results as a function of h~ 2 . The one-mode 
high-field limits i? Cj00 and fc CjOC are shown in the figures 
as open triangles. The corresponding numerical many- 
mode results are given as solid triangles. The experimen- 
tal data for R c and k c are consistent with the expected 
dependence on h, but respectively fall about 4% above 
and 1% below the calculation. The data for R c and k c at 
the highest experimental field h ~ 80 are already within 
about 6% and 1% respectively of the infinite-field value. 
Thus it seems unlikely that qualitatively new phenom- 
ena could be discovered by measurements at even higher 
fields. 
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FIG. 22. The critical wavenumber as a function of 1/h 2 . 
The solid line was calculated using the parameters of Ref. 

The open triangle is k Ct00 from Eq. [l3|. The solid trian- 
gle is the many-mode numerical result for k Ct00 . The vertical 
bars indicate the location of the codimension-two point (right 
bar) and the tricritical point (left bar). The dashed line shows 
what happens to the theoretical curve if the viscosity 04 is 
increased by 7.5% (the data for k c remain unchanged). 

We examined whether the small difference between the 
theory and the experiment could be removed by small 
adjustments in the values of the fluid properties. We 
found that an increase by 7.5% of 014 yielded the plusses 
and crosses for the data in Fig. |2l] and the dashed lines 
in Figs. [2l] and [2^ (the data for k c in Fig. [2^ are not 
affected by changing 04). The adjustment of 04 pro- 
duced an excellent fit for both k c and R c . However, it 
spoiled the excellent agreement for R c along the oscilla- 
tory branch below the codimension-two point shown in 
Fig. H and did not significantly reduce the difference be- 
tween calculation and experiment for k c at small h which 
is shown in Fig. ||. Various other attempts to adjust the 
fluid properties used in the theoretical calculations were 
unsuccessful in yielding improved overall agreement be- 
tween theory and experiment over the entire field range. 



F. Nonlinear States at High Fields 



Beyond the codimension-two point the finite- 
amplitude flow was split into three regions distinguished 
by their characteristic wavenumbers k n ~ 3.4, 3.6 and 
3.9. These regions are shown in Fig. |l5|. Figure ^ shows 
some characteristic patterns. Qualitatively the patterns 
appear similar, each exhibiting spatio-temporal chaos. 
The transitions between these state depended on both h 
and R. They were determined by measuring the change 
in k n as h was varied at a given Rayleigh number. As 
already mentioned above, to our knowledge there are no 
theoretical predictions for these patterns. 
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FIG. 23. Time sequences of the spatially and temporally 
chaotic flow for R — 4000 and at different values of h for cell 
5. The images in each row were taken at the field indicated in 
the leftmost image. Time increases from left to right. Images 
were taken in one hour intervals. The wavenumbers of the 
patterns are h = 47: k n = 3.4, h — 56: k n = 3.6, h = 57.5: 
k n = 3.9. 
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